UTHEP-348 
UTCCP-P-19 
February, 1997 

Scaling in SU(3) Pure Gauge Theory 

with a Renormalization Group 
Improved Action 

Y. Iwasakif' b K. Kanayaf' b T. Kanekof and T. Yoshie a ' b 

a Institute of Physics, University of Tsukuba, Ibaraki 305, Japan 
b Center for Computational Physics, University of Tsukuba, 
Ibaraki 305, Japan 

We study the scaling properties of the static quark potential and the ratio 
of the critical temperature T c to the square root of the string tension a in the 
SU(3) pure gauge theory using a renormalization group improved action. We 
first determine the critical coupling (3 C on lattices with temporal extension 
N t = 3, 4, and 6, and then calculate the static quark potential at the critical 
couplings on lattices at zero temperature. We note that the static quark 
potentials obtained are rotationally invariant with errors of at most 1 - 2 % 
in all the three cases, and that the potential V(R) in physical units scales in 
the whole region of R investigated. The values of T c j yfa for the three cases 
in the infinite volume limit are identical within errors. We estimate the value 
in the continuum limit to be T c j ' yfo = 0.656(4), which is slightly larger than 
the value in the continuum limit from the one-plaquette action, 0.629(3). 



1 Introduction 



In numerical studies of lattice QCD, it is important to control and reduce 
finite lattice spacing effects. Several improved actions have been proposed 
for this purpose and some of them have been tested for the scaling behavior 
of the critical temperature T c of the finite temperature deconfining transition 

In this work we study the scaling properties of the static quark potential 
and the ratio of the critical temperature to the square root of the string 
tension a, T c /y/a, in the SU(3) pure gauge theory, using a renormalization 
group (RG) improved action |J: 

S T = f ( c o EC 1 >< 1 lo °P) + c i EC 1 >< 2 lo °P)} (!) 

with c\ = —0.331 and Co = 1 — 8c\, where (3 = 6/g 2 (g is the gauge cou- 
pling). In Eq.(|lD, the loops are defined by the trace of ordered product of 
link variables and each oriented loop appears once in the sum. 

This paper is organized as follows. First we determine the critical coupling 
/3 c 's for the finite temperature deconfining phase transition on 9 3 x 3, 12 3 x 4 
and 18 3 x 6 lattices in Sec. 2. We also perform simulations on 12 3 x 3, 15 3 x 3, 
16 3 x 4, and 20 3 x 4 lattices for a finite size scaling study. Then the quark 
potentials at the three /3 c 's are calculated from smeared Wilson loops on 
9 3 x 18, 12 3 x 24 and 18 3 x 36 lattices, respectively, in Sec. 3. The string 
tension is extracted from the quark potential assuming that the potential 
takes a form of a sum of a Coulomb term and a linearly rising potential. In 
Sec. 4, scaling behavior of the quark potential and that of the ratio T c /y/a 
are examined. Finally, the value of the ratio T c j \fa in the continuum limit 
and in the infinite volume limit is estimated. 

2 Critical coupling j3 c 

In order to determine the critical coupling /3 C for the finite temperature phase 
transition, we perform simulations on 9 3 x 3, 12 3 x 4 and 18 3 x 6 lattices. 
The critical temperature T c is given by T c = l/(aN t ), where N t is the linear 
extension of the lattice in the temporal direction and a is the lattice spacing 
at the critical coupling. Note that the physical spatial volumes are identical 



1 



for all the three cases, (N s a) 3 = (3N t a) 3 = (3/T c ) 3 , where N s is the linear 
extension of the lattice in the spatial direction. 

We also perform simulations on lattices with different spatial volumes 
for an estimation of the infinite volume limit of (3 C using finite size scaling 
analyses. The previous results for the case of the standard one-plaquette 
action on spatially large lattices |7|, |8|] indicate that extrapolations from small 
lattices with the aspect ratio N s /N t < 3 result in sizable systematic errors in 
the values of (3 C in the infinite volume limit. Therefore, we restrict ourselves 
to lattices N s /N t > 3 in this paper. We perform simulations on 12 3 x 3, 
15 3 x 3, 16 3 x 4, and 20 3 x 4 lattices for finite size analyses. We reserve finite 
size study of N t = 6 lattices for future investigation. 

Gauge fields are updated by a Cabibbo-Marinari-Okawa pseudo heat bath 
algorithm with 8 hits both for the simulations at finite temperatures and at 
zero temperature discussed in the next section. The simulation parameters 
are compiled in Table |1[ We measure Wilson loops and Polyakov line every 
10 sweeps. Their expectation values are summarized in Tables |] - |[ (For 
the deconfinement fraction, see below.) 

The values of the critical coupling f} c are determined as the peak location 
of the susceptibility x °f the Z(3) rotated Polyakov line Q: 



The results of the susceptibility calculated using the spectral density 
method § on the 9 3 x 3, 12 3 x 4, and 18 3 x 6 lattices are shown in Fig. |T[ 
The results obtained at several simulation points are consistent with each 
other within the errors and form a clear peak structure. The value of /3 C is 
determined from the data at the f3 which is the closest to f3 c . The errors are 
estimated using a single-elimination jack knife method. The bin size in the 
jack knife method is determined by investigating the bin size dependence of 



(2) 




(3) 



(4) 
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the errors of Q, shown in Fig. 0. We note that the jack-knife errors of /3 c 's 
are stable for the bin size larger than those adopted, as shown in Fig. ^. The 
values of /3 c 's and their jack-knife errors are summarized in Table []. 

There are several alternative definitions of (3 C on finite lattices. A popular 
method is to measure the "deconfinement fraction" r given by r = (3p — 1)/2 
where p is the probability such that arg P £ (— 7r/9, 7r/9), (27r/3 — 7r/9, 27r/3 + 
7r/9), or (— 27r/3 — 7r/9, — 27r/3 + 7r/9), and to define (3 C as a point where r 
takes a given value. Our results of r as a function of f3 for the case of the 
aspect ratio N s /N t = 3 are shown in Fig. [|. See also Tables |]-|. We find that 
the deconfinement fraction r is approximately 0.75 at f3 c determined from the 
peak location of the susceptibility, as summarized in Table |9|. We note that 
this fact for the deconfinement fraction is also realized in data || obtained for 
the standard one-plaquette action on large lattices with high statistics (see 
Table ITO). The condition r = 3/4 is the criterion taken in Ref. [ITU] for the 



determination of (3 C . (See also the discussions in Refs. [[TTJ, |I2"]|.) However, the 
volume dependence of the corrections of f3 c to the infinite volume limit is not 
known in this case.Q On the other hand, a scaling relation is well established 
for the 13 C determined from the peak location of the susceptibility. Therefore, 
we concentrate on (5 C determined from the peak location of the susceptibility 
for finite size scaling analyses. 

In the following, we denote the j3 c on the 9 3 x 3, 12 3 x 4, and 18 3 x 6 
lattices as /3 C (9 3 x 3), /3 C (12 3 x 4), and /? c (18 3 x 6), respectively. 



1 The value r — 3/4 corresponds to the case that the four peaks of the histogram 
of P in the complex plane have the same volume fraction [fl2|| , assuming uniformity of 
the distribution in terms of argP in the confining phase. For the g-state Potts models 
with large q, the value of (3 C which corresponds to the case where {q + 1) peaks have 
the same volume fraction is shown to yield the correct infinite volume value of /3 C up to 
exponentially suppressed corrections |L3|. However, in the SU(3) gauge theory, uniformity 
of distribution in terms of argP in the confining phase is not well satisfied. Therefore, 
r = 3/4 does not strictly correspond to the case of equal weight of four peaks. Thus, in 
contrast to the case of (3 C from the peak location of the susceptibility, no rigorous scaling 
relation is known for the (3 C determined from the deconfinement fraction. In practice, when 
we adopt (3 C determined from r = 3/4 and assume either a linear volume dependence or 
an exponential volume dependence, we obtain a result for T c j yfa in the continuum limit 
which agrees, within errors, with that derived in the text using /3 C from the peak location 
of the susceptibility. 



3 



3 String tension 



We evaluate the string tensions at f3 = f3 c (9 3 x 3), /? c (12 3 x 4), and /5 C (18 3 x 6) 
on lattices at zero temperature: 9 3 x 18, 12 3 x 24, and 18 3 x 36 lattices, 
respectively. Note that the spatial sizes of the lattice are the same as those 
for the finite temperature simulations in all the three cases. The ratio N t /N s 
is also fixed to 2. The simulation parameters are summarized in Table |TT[ 
After thermalization sweeps, we measure Wilson loops every 200 sweeps. The 
spatial paths of the loops are formed by connecting one of the spatial vectors 
shown in Fig. ||. 

In order to extract the ground state contribution to the potential, we 



adopt the smearing technique proposed in Ref . ||14|| : Each spatial link U is 
replaced with an SU(3) matrix U ncw which maximizes ReTr [FU new ], with F 
being the sum of the spatial staple products of link variables around U. We 
perform this procedure up to 10, 30 and 40 steps on the 9 3 x 18, 12 3 x 24 and 
18 3 x 36 lattices, respectively. Measurements are carried out every smearing 
step on the 9 3 x 18 and every 2 smearing steps on the other lattices. With 
this smoothing procedure the behavior of the effective mass 

m eff = log [W(R, T)/W(R, T + 1)] (5) 

in terms of T is much improved, especially for large R as shown in Fig. || 

In the following, we discuss separately the results of the potential V(R) 
at /3 C for N t = 4 and 6, and that for N t = 3, because in the former case we 
are able to extract the coefficient of the Coulomb term by a straightforward 
fitting procedure with examining the stability of the fit, while in the latter 
case it is hard to fix it solely from the data due to a small number of the data 
points caused by the coarseness of the lattice at (3 C (9 3 x 3) (see discussions 
below). 



3.1 Results at /3 C (12 3 x 4) and /3 C (18 3 x 6) 

The potential V(R) and the overlap function C(R) are extracted by a fully 
correlated fit of Wilson loops to the form 

W(R, T) = C(R)- exp [-V(R)-T] . (6) 

The fitting range is determined by examining carefully x 2 / df and stability 
of V(R) against the smearing step. Fig. |7| shows the results of x 2 /df and 
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V(R) versus the smearing step at R = 4.0 for the case of /3 C (12 3 x 4). When 
we take the fitting range T = 3 — 5, we find that x 2 /df < 1 and V(R) is 
quite stable after four smearing steps, while the choice of the fitting range 
T = 2 — 4 leads to X 2 /df much larger than 1 and a significant variation 
of V(R) against the smearing step. We find that the choice of the fitting 
range T = 3 — 5 leads to reasonable x 2 /df and stability of V(R) against 
the smearing step for all R except 2 a/6 (where x 2 /df takes a little large 
value ~ 2.5, though the stability is satisfied). This stability implies that the 
contamination from excited states is negligible small. Therefore, we take the 
fitting range T = 3 — 5 for the data at /3 C (12 3 x 4). The T range 4 — 7 at 
/? c (18 3 x 6) is determined in a similar way. 

We determine the optimum number of smearing steps for each R in such 
a way that C(R) takes the largest value under the condition C(R) < 1 which 
we call the "optimum smearing step". We note that x 2 /df is stable (< 
1) against a variation of the smearing step when C(R) ~ 1. The optimum 
smearing steps thus determined are about 8 at /? c (12 3 x4), and are distributed 
from 12 to 40 at /3 C (18 3 x 6) (see Tables [12] and [T|).[] We take the value of 
V(R) at the optimum smearing step. The systematic error due to the choice 
of the smearing step is much smaller than the statistical error, because the 
value of V(R) is stable against the smearing step as mentioned above, and 
therefore we neglect it in the following. 



The values for V(R) are summarized in Tables [12] and 13. Statistical 



errors are estimated by the jack-knife method with bin size 1. Note that 
measurements are performed every 200 sweeps. We confirm that the errors 
are quite stable against the bin size. 

The string tension is determined by fitting V(R) to the rotationally in- 
variant ansatz 

rv 

V(R) =V - — + a lat R, (7) 
ri 

where <j\ at = &a 2 is the string tension in lattice units. We take into account 
the correlations among V(R) at different R using the error matrix derived 

2 Wc find that the value of C(R) for R = 1.0 - 2.0 on the 18 3 x 36 lattice is greater 
than 1 at all smearing steps < 40. We have checked using 20 configurations that more 
than 60 smearing steps are needed to get C{R) < 1 for these i?'s. Because we do not use 
these small loops for the fit of the potential, we stop the smearing steps at 40 times. 
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from those for W(R, T). The fitting ranges we take are 



R=V6-AV2 (12 3 x24), 

R=2y/3 - 4v/5 (18 3 x36). 1 ' 

These ranges (-R mm - -Rmax) are determined by investigating the stability of 
fits and the value of x 2 /df as explained in the following. As we increase -R mm , 
instability of the fit first appears in the result of a, while the results of V and 
0i at are stable. The error of a becomes abruptly large as R m i n increases: e.g. 
at /3 C (12 3 x 4) with ,R max = 4^ fixed, a = 0.332(11), 0.295(14), 0.154(101), 
and -0.040(121) for -R m i n = a/5, \/6, 2\/2, and 3.0, respectively. Therefore, we 
restrict candidates for R m i n to those for which the error of a is less than 50% 
of the central value. We find that x 2 /df is stable and ~ 1 for 5 < R max < 6 
[7 < i? max < 9 ] at /3 C (12 3 x 4) [/? c (18 3 x 6)] which we take as the candidates 
for i? max . The fitting range is determined by the condition that x 2 /df takes 
a value nearest to 1 in all the combinations of the candidates for -R m i n and 
R max . The values of x 2 /df are 1.5 and 1.2 at /5 C (12 3 x 4) and /3 C (18 3 x 6), 
respectively, for the -R m i n and -R max adopted. We have checked that the results 
of a and o"i at are stable for all candidates of (-R m i n , -Rmax) which satisfy 

2.0 < R min < y/E, 5.0 < i? max < 6.0 (12 3 x 24) 

3.0 < R min < 3V2, 7.0 < i? max < 9.0 (18 3 x 36). 1 ' 

Note that the changes of the fitting ranges of R at these two /3's are con- 
sistent with the change of the scale between (3 = /3 c (12 3 x4) and /3 c (18 3 x6), 
that is, the ratio of 4 to 6. 

The results of Vo, a, 0"i at , and their jack knife errors are summarized in 
Table 0. The values of V(R) are plotted in Fig. [8], where different symbols 
correspond to different units of spatial path of Wilson loops. The values of 
V(R) obtained from six types of Wilson loops are excellently fitted to the 
rotationally invariant form, Eq.(|7|). The deviations of the data at /9 C (12 3 x 4) 
from the fitted curve are less than 2% and the average of them is about 0.4%. 
For the data at /3 C (18 3 x 6), the deviations are at most 1% with the average 
about 0.3%. 

We note that the results of a are consistent with a constant within the 
errors. The resulting a ~ 0.296 is slightly larger than 7r/12 ~ 0.262 derived 



in a string model |T5[]. We also perform fits with the value of a fixed to 



7r/12. Then the values obtained are cri at = 0.1527(14) and 0.0667(6) at 
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/? c (12 3 x 4) and /3 C (18 3 x 6), respectively. The values for the ratio T c /^/a 
using these results are consistent with our final results using the values in 
Table IT4] within one standard deviation. 



3.2 Results at &(9 3 x 3) 

We obtain the potential V(R) at (3 C (9 3 x 3) by fitting W(R,T) to the form 
Eq.(|6|) with the fitting range T = 2 — 4 . The fits with this fitting range 
have desirable properties similar to those at the other two /3 c 's discussed in 
the preceding subsection; reasonable x 2 /df and stability of V(R) against the 
smearing step. 

When we make a fit of the potential to the form Eq. (|7]) , we find that the 
-Rmin dependence of a is stronger than the cases discussed in the previous 
subsection, while the fits are quite stable against -R max like in the previous 
cases. This is due to the fact that we have only small number of data points 
at small R caused by the coarseness of the lattices at /3 C (9 3 x 3). Therefore a 
small deviation from the rotational invariance at R = R m i n sometimes affects 
the value of a sizable. As a result, we are not able to find an -R m i n region for 
which a is stable. 

Therefore, we perform two kinds of fits at (3 — (3 C (9 3 x 3): In the first 
fit, we fix the value of a to the average value 0.296 of those at the other 
two /3 c 's which are constant within the errors. We set the fit range to be 
R = 2 — 2v^5 so that the physical R range is consistent with the ranges 
at /3 c (12 3 x4) and /3 c (18 3 x6). As shown in Fig. ||, the fit well reproduces the 
data even at R < R m m- m the other fit, we perform fit without fixing the 
value of a, for the ranges -R m i n = V2, V3, and 2 and -R max = 2y/E. These 
values of -R m i n in physical units correspond to those at the other two /3 c 's for 
which the stability of a is observed. 

We take the results of the former fit with a fixed as the central values of 
o"iat and Vq. The statistical errors are obtained by the jackknife method with 
bin size 1. We then take the upper bounds and lower bounds of Vq and cr lat 
obtained by the fits a unfixed, as systematic errors. The results of Vq and 
cxiat with the errors are given in Table 0. The potential data are shown in 
Fig. ^ together with its fit curve (a fixed to 0.296). The deviations from the 
fit are at most 2% and the average of them is about 0.5%, which indicates 
that the rotational invariance is well restored even at this small value of (3. 

We also perform a fit with a fixed to n/12 to find o"i at = 0.2607(16). 
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The ratio T c j yfo using this result is consistent with our final result using the 
value in Table [IB] within the errors. 



4 Scaling properties 

In Fig. [LO], the values of V v ]y YS /T c are shown as a function of Ra ■ T c , where 
Vphys = V/a is the potential in physical units. We note that the data on all 
the lattices are in excellent agreement in the whole Ra-T c region. This implies 
scaling of our potential data in the range of (3 values investigated. It might 
be emphasized again that the deviation of the data from the rotationally 
invariant fit is at most 2 % for the N t = 3 and 4 cases and 1 % for the N t = 6 
case. 

Using the results presented in the preceding section, we obtain the values 
of T c j ^fo on the lattices with finite spatial volume 9 3 , 12 3 , and 18 3 , which 
equal to (3/T c ) 3 m (2.2fm) 3 in physical units: 

f 0.660(3)0 (iV* = 3) 
T c /v^(finite volume) = < 0.647(5) (N t = 4) (10) 

[ 0.651(6) {N t = 6). 

The number in the first brackets is the statistical error and the second one 
for Nt = 3 is the systematic error due to uncertainty of the fitting range. 

In order to estimate the values of T c f \fa in the infinite volume limit, we 
first obtain finite size scaling relations |7], |8j 

P c (N t , iV 3 ) = p c (N t , oo) - 0.122(54) • N t 3 /N 3 s (N t = 3), (11) 

and 

P e (N u iV 3 ) = p c (N t , oo) - 0.133(63) ■ N t 3 /N* (N t = 4), (12) 

from the data of /3 C on the N s /N t = 3, 4, and 5 lattices (see Fig. [TTJ) . We note 
that the slopes of /3 c (N t , iV 3 ) in Nf /N 3 in the two relations are independent 
of N t within the errors, as observed previously in the case of the standard 
one-plaquette action ||. Therefore, we assume the relation ( |I~2"D also for 
N t = 6. Then we have 

f 2.1551(12) (N t = 3) 
f3 c (N t , oo) = \ 2.2879(11) {N t = 4) (13) 
[ 2.5206(30) {N t = 6). 
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The values of the string tension at (3 c {N t , oo) are estimated assuming an 
exponential scaling of y/d^t in terms of j3 [|l^]. We obtain 

= 25.26(43) ■ exp [ -1.828(7)-/? ] , (14) 

by fitting the values of ai at at /3 C (9 3 x 3), /3 C (12 3 x 4), and /3 C (18 3 x 6) as 
shown in Fig. [IJ| This relation is used to compute the shifts in oi at from the 
values at f3 c (N t , N 3 ) to those at /3 c (N t , oo). The values of o"i at (/3 c (A^, oo)) are 
obtained by adding the shifts to those of eri a t(/3 c (iV t , Nf)) given in Tables [14| 
and [TBI: 



0.2516(26)(11)(±I? ) (N t = 3) 
a lab (p c (N t ,oo)) = { 0.1464(25)(6) (N t = 4) (15) 

0.0644(12)(7) {N t = 6). 

The number in the first brackets is the statistical error, the second one is 
the error due to the error in the values of /3 c (N t , oo), and the third one for 
N t = 3 is the systematic error due to uncertainty of the fitting range. 
Finally, we obtain 

f 0.665(3X1)0 (N t = 3) 
T c /v^(infinite volume) = < 0.653(6) (1) (N t = 4) (16) 

[ 0.657(6)(4) (N t = 6). 

The origins of the errors are the same as in Eq. flToT) . Our three values are 
consistent with a constant within the errors. A weighted average of the values 
given in Eq. fllB]) gives 

T c /y/a = 0.656(4) (17) 



in the continuum limit. 

Using the experimental value a = <Ji a t/a 2 = (420MeV) 2 , we obtain a ~ 
0.23, 0.18, and 0.12 fm at f3 c for Nt = 3, 4, and 6, respectively. Thus 
the scaling behavior for the ratio T c j \J~d starts at least around a ~ 0.23fm 
with the RG improved gauge action. From Eq.(|T7D we also obtain T c pa 
276(2)MeV. 



Our results ( |T6D are shown in Fig. |I3| together with the results using 
other actions j|, 0. Our result T c j ' yfo = 0.656(4) in the continuum limit is 
slightly larger than the value with the standard action 0.629(3) [16|]. We also 
compare our results with those derived from the torelon mass fi{L) which 
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is calculated from Polyakov line correlators on a lattice of spatial size L. 



Defining u(L) = fi(L)/L, we extrapolate the values of T c /ya(L) to the 
continuum limit. Then the value of T c / ■sfo is estimated assuming the relation 
a = (x(L)+7r/ (3L 2 ) derived in a string model JL7]]. (We neglect the corrections 
due to the shift f3 c (V = oo) — f3 c {V = N%).) For a fixed point action 0, we 
obtain T c /y/a = 0.617(5) using the data for N t = 2, 3 and 4 with N s = 2N t . 
The result is about 6% smaller than our result fllTI). For a tadpole-improved 
Symanzik action ||, we obtain T c /yfa = 0.649(5) using the data for N t = 3 
and 4 with N s = 2N t . The result is consistent with our result. 

Numerical simulations are performed with Fujitsu VPP500/30 and HI- 
TAC H6080-FP12 at the University of Tsukuba. We thank Akira Ukawa for 
valuable discussions. This work is in part supported by the Grants-in-Aid of 
Ministry of Education, Science and Culture (Nos.07NP0401, 07640375 and 
07640376) and the University of Tsukuba Project Research in 1996. 
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Table 1: Parameters of finite temperature simulations 



1 J-X * 

lattice size 


(3 


sweep 


therm. 


9 3 x 3 


2.125 


100000 


30000 




O 1 A C 

Z.145 


100000 


oOOOO 




2.160 


70000 


30000 


12 3 x 3 


2.150 


100000 


50000 




2.155 


100000 


40000 


15 3 x 3 


2.150 


180000 


80000 


12 3 x 4 


2.250 


12000 


2000 




2.275 


125000 


40000 




2.300 


10000 


1500 


16 3 x 4 


2.283 


220000 


40000 




2.290 


240000 


40000 


20 3 x 4 


2.2875 


270000 


80000 


18 3 x 6 


2.5000 


120000 


15000 




2.5125 


256000 


50000 




2.5250 


210000 


60000 




2.5375 


135000 


5000 


Table 2: Results obtained on the 9 3 x 3 lattice. 




13 = 2.125 


13 = 2.145 


13 = 2.160 


lxl Wilson loop 


0.575350(61) 


0.58207(12) 


0.58752(15) 


1x2 Wilson loop 


0.32181(11) 


0.33120(22) 


0.33917(28) 


2x2 Wilson loop 


0.10721(18) 


0.11600(32) 


0.12436(43) 


Polyakov line 


0.0675(23) 


0.1055(40) 


0.1575(52) 


deconfinement fraction 


0.448(23) 


0.666(29) 


0.895(29) 
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Table 3: Results obtained on the 12 3 x 3 lattice. 





P = 2.150 


P = 2.155 


lxl Wilson loop 


0.58329(16) 


0.58546(11) 


1x2 Wilson loop 


0.33277(31) 


0.33607(20) 


2x2 Wilson loop 


0.11716(47) 


0.12094(31) 


Polyakov line 


0.0972(74) 


0.1300(38) 


deconfinement fraction 


0.691(53) 


0.899(24) 



Table 4: Results obtained on the 15 3 x 3 lattice. 



P = 2.150 

0.58321(17) 

0.33260(32) 

0.11688(48) 

0.0862(57) 

0.715(52) 



lxl Wilson loop 
1x2 Wilson loop 
2x2 Wilson loop 
Polyakov line 
deconfinement fraction 



Table 5: Results obtained on the 12 3 x 4 lattice. 





P = 2.250 


P = 2.275 


P = 2.300 


lxl Wilson loop 


0.608085(91) 


0.614037(60) 


0.62007(28) 


1x2 Wilson loop 


0.36552(18) 


0.37384(12) 


0.38255(52) 


2x2 Wilson loop 


0.14448(25) 


0.15257(19) 


0.16175(76) 


Polyakov line 


0.0374(25) 


0.0651(31) 


0.1213(16) 


deconfinement fraction 


0.342(33) 


0.640(23) 


0.981(23) 
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Table 6: Results obtained on the 16 3 x 4 lattice. 





P = 2.283 


P = 2.290 


lxl Wilson loop 


0.615677(33) 


0.617594(35) 


1x2 Wilson loop 


0.376040(65) 


0.378940(69) 


2x2 Wilson loop 


0.15451(10) 


0.15781(11)) 


Polyakov line 


0.0549(30) 


0.0906(30) 


deconfinement fraction 


0.583(39) 


0.853(27) 



Table 7: Results obtained on the 20 3 x 4 lattice. 



P = 2.2875 

0.616941(56) 

0.37794(11) 

0.15663(19) 

0.0744(57) 

0.768(48) 



lxl Wilson loop 
1x2 Wilson loop 
2x2 Wilson loop 
Polyakov line 
deconfinement fraction 



Table 8: Results obtained on the 18 3 x 6 lattice. 



P = 2.5000 /3 = 2.5125 p = 2.5250 /3 = 2.5375 



lxl Wilson loop 
1x2 Wilson loop 
1x3 Wilson loop 
2x2 Wilson loop 
2x3 Wilson loop 
3x3 Wilson loop 
Polyakov line 
deconfinement fraction 



0.655687(59) 

0.431669(40) 

0.288927(53) 

0.208714(77) 

0.109884(87) 

0.049770(87) 

0.0328(32) 

0.555(51) 



0.657691(11) 

0.434539(25) 

0.291912(33) 

0.211717(46) 

0.112330(53) 

0.051480(54) 

0.0409(21) 

0.645(36) 



0.659676(11) 

0.437385(24) 

0.294884(33) 

0.214735(47) 

0.114813(54) 

0.053261(58) 

0.0559(22) 

0.861(31) 



0.661649(12) 

0.440230(24) 

0.297874(32) 

0.217780(34) 

0.117350(55) 

0.055096(59) 

0.0691(21) 

0.960(15) 
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Table 9: The value of (3 C and its jack knife error determined by the peak 
location of the susceptibility, and the deconfmement fraction r at (3 C , together 
with the bin size adopted. 





Pc 


r 


bin size 


9 3 x 3 


2.1508(12) 


0.757(25) 


1000 


12 3 x 3 


2.1528(9) 


0.771(48) 


3000 


15 3 x 3 


2.1546(11) 


0.894(33) 


8000 


12 3 x 4 


2.2827(16) 


0.774(47) 


3000 


16 3 x 4 


2.2863(10) 


0.765(37) 


6000 


20 3 x 4 


2.2865(9) 


0.742(52) 


10000 


18 3 x 6 


2.5157(7) 


0.698(34) 


3000 



Table 10: Deconfinement fraction r at (3 C determined from the susceptibility 
in the case of the standard one-plaquette action. Data of the Polyakov line 
are taken from a previous simulation by the QCDPAX Collaboration ||. 



lattice 


Pc 


r 


12 3 x 24 x 4 


5.69149 


0.790(12) 


24 3 x 36 x 4 


5.69245 


0.732(46) 


20 3 x 6 


5.8924 


0.805(26) 


24 3 x 6 


5.89292 


0.786(27) 


36 3 x 48 x 6 


5.89379 


0.739(40) 
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Table 11: Simulation parameters for the quark potential measurements. 



lattice size 


P 


thermalization 


# of conf. 


9 3 x 18 


2.1508 


5000 


400 


12 3 x 24 


2.2827 


5000 


200 


18 3 x 36 


2.5157 


10000 


100 



Table 12: Potential V(R) and overlap function C(R) at /9 c (12 3 x4) obtained 
on the 12 3 x 24 lattice. The optimum smearing step N opt is also given. 



R 


unit of R 


V(R) 


C{R) 


N opt 


1.000 


(1,0,0) 


0.47408(35) 


0.9997(8) 


10 


1.414 


(1,1,0) 


0.63207(60) 


0.9987(14) 


10 


1.732 


(1,1,1) 


0.72136(89) 


0.9955(23) 


10 


2.000 


(1,0,0) 


0.7811(12) 


0.9916(31) 


12 


2.236 


(2,1,0) 


0.8304(11) 


0.9880(25) 


10 


2.449 


(2,1,1) 


0.8759(12) 


0.9977(29) 


8 


2.828 


(1,1,0) 


0.9472(17) 


0.9990(43) 


8 


3.000 


(1,0,0) 


0.9843(27) 


0.9857(72) 


10 


3.000 


(2,2,1) 


0.9806(19) 


0.9799(48) 


8 


3.464 


(1,1,1) 


1.0574(32) 


0.9739(86) 


8 


4.000 


(1,0,0) 


1.1507(15) 


0.9779(41) 


8 


4.243 


(1,1,0) 


1.1926(45) 


0.982(12) 


8 


4.472 


(2,1,0) 


1.2317(41) 


0.976(11) 


8 


4.899 


(2,1,1) 


1.3053(55) 


0.971(15) 


8 


5.000 


(1,0,0) 


1.3222(22) 


0.9881(63) 


8 


5.196 


(1,1,1) 


1.34430(55) 


0.9925(57) 


6 


5.657 


(1,1,0) 


1.426(11) 


0.982(29) 


8 


6.000 


(1,0,0) 


1.479(19) 


0.978(53) 


8 


6.000 


(2,2,1) 


1.49012(78) 


0.9966(66) 


6 
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Table 13: The same as Table [TJ at /3 c (18 3 x6) obtained on the 18 3 x 36 lattice. 
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V{H) 
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O.UUU 
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k i nf; 
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K RK7 
0.00 ( 


f i 1 n"\ 
^1,1, U J 


U.y40Z^Zo ) 


n ofiO/i ^07^ 
u.ysy4^y / j 


1 Q 

lo 


6.000 


(1,0,0) 


0.9644(39) 


0.993(15) 


14 


6.000 


(2,2,1) 


0.9721(29) 


0.992(10) 


16 


6.708 


(2,1,0) 


1.0254(34) 


0.982(11) 


18 


6.928 


(1,1,1) 


1.0374(43) 


0.984(15) 


16 


7.000 


(1,0,0) 


1.04334(53) 


0.9905(64) 


16 


7.071 


(1,1,0) 


1.0501(18) 


0.9983(84) 


16 


7.348 


(2,1,1) 


1.0674(41) 


0.995(14) 


14 


8.000 


(1,0,0) 


1.11216(68) 


0.9953(78) 


14 


8.485 


(1,1,0) 


1.15463(52) 


0.9998(84) 


16 


8.660 


(1,1,1) 


1.1714(32) 


0.988(14) 


18 


8.944 


(2,1,0) 


1.1743(22) 


0.9870(78) 


12 


9.000 


(1,0,0) 


1.17869(54) 


0.9794(86) 


14 


9.000 


(2,2,1) 


1.18690(47) 


0.9975(78) 


12 
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Table 14: Results of Vo, a, and oi a t at /3 c (12 3 x4) and /3 c (18 3 x6) obtained on 
the 12 3 x24 and 18 3 x36 lattices, respectively. 









V 


a 


C"lat 


/3 C (12 3 


x 4) 




0.630(20) 


0.295(14) 


0.1493(25) 


/3 C (18 3 


x 6) 




0.627(18) 


0.297(19) 


0.0655(12) 


Table 15: 


rhe same as 


Table 12 at [3 C (9 3 x 


3) obtained on the 9 3 


x 18 lattice. 


R 




unit of R V(R) 


C(R) 


Nopt 


1.000 




(1,0,0) 


0.54892(35) 


0.9982(5) 


4 


1.414 




(1,1,0) 


0.75390(63) 


0.9984(10) 


4 


1.732 




(1,1,1) 


0.87752(99) 


0.9942(16) 


4 


2.000 




(1,0,0) 


0.9601(13) 


0.9981(22) 


5 


2.236 




(2,1,0) 


1.0349(13) 


0.9965(21) 


4 


2.449 




(2,1,1) 


1.1049(15) 


0.9963(25) 


2 


2.828 




(1,1,0) 


1.2101(22) 


0.9945(39) 


4 


3.000 




(1,0,0) 


1.2640(31) 


0.9836(55) 


5 


3.000 




(2,2,1) 


1.2685(23) 


0.9591(39) 


3 


3.464 




(1,1,1) 


1.3890(41) 


0.9736(73) 


4 


4.000 




(1,0,0) 


1.5420(18) 


0.9750(34) 


5 


4.243 




(1,1,0) 


1.6001(60) 


0.987(12) 


4 


4.472 




(2,1,0) 


1.6634(62) 


0.978(11) 


4 


4.899 




(2,1,1) 


1.7785(53) 


0.973(13) 
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Table 16: Results of Vo and cri at at /3 C (9 3 x 3) obtained on the 9 3 x 18 lattice. 
The numbers in the first brackets are statistical errors and the second are 
systematic errors due to uncertainty of the fitting range. 



Vo 


Clat 


0.598(60)(i^) 


0.2554(26) (i^ ) 
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Figure 1: Susceptibility x of the Z(3) rotated Polyakov line Q on the 9 3 x3, 
12 3 x4 and 18 3 x6 lattices. Each curve is obtained from the data at the 
simulation point indicated by the filled symbol using the spectral density 
method. Errors for open symbols are computed by a jack-knife method 
combined with the spectral density method. 
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Figure 2: Bin size dependence of the jack knife error of f2. 
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Figure 3: Bin size dependence of the jack knife error of f3 c 
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Figure 4: Deconfiriement fraction r on the 9 3 x3, 12 3 x4 and 18 3 x6 lattices. 
Each curve is obtained from the data at the simulation point indicated by 
filled triangles. The values of (3 C determined from the peak position of the 
susceptibility are shown by filled circles. 
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Figure 4: continued 
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Figure 5: Units of spatial paths of Wilson loops. 
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Figure 6: T dependence of m cS = log [W(R, T) /W(R, T + 1)] without (a) 
and with (b) smearing on the 12 3 x 24 lattice. The results shown in (b) are 
obtained with the optimum smearing step explained in the text. 
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Figure 7: Smearing step dependence of x 2 /df and V(i2) at R = 4.0 on the 
12 3 x 24 lattice. 
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Figure 8: Data for the potential and its fitting curves on the 12 3 x24 and 
18 3 x36 lattices. The values of a are determined from a = (420MeV) 2 . The 
legends for the symbols represent the units of spatial paths of Wilson loops 
shown in Fig. ||. 
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Figure 10: Vphys/^c vs. Ra-T c . The constant term in the potential is fixed so 
that the potentials have the same value at Ra-T c = 1.0. 
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Figure 11: Finite size scaling of f3 c on N t = 3 and 4 lattices. The solid lines 
are the results of a linear fit and the dashed lines indicate the location of the 
infinite volume limit. 
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Figure 12: Square root of the string tension v^iat a ^ A= on the finite volume 
lattices simulated together with its fitting curve. 



31 



0.750 




• RG improved action 






o standard action 






a (1 ,2) Symanzik action 






[> (2,2) Symanzik action 


0.700 




V (1 ,2) Symanzik action (tadpole improved) 


0.650 


« 




















0.600 




5 ; 





0.00 0.05 0.10 

1/Nt 2 



Figure 13: T c / 'y/a vs. l/N t 2 . All symbols represent the values in the infinite 
volume limit. The errors for our data in the cases of N t = 4 and 6 are statis- 
tical, while those for N t = 3 include the systematic error due to uncertainty 
of the fitting range. 
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